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Abstract. Fractals define a new and interesting realm for a discussion of ba- 
sic phenomena in quantum field theory and statistical mechanics. This interest 
results from specific properties of fractals, e.g., their dilatation symmetry and 
the corresponding absence of Fourier mode decomposition. Moreover, the ex- 
istence of a set of distinct dimensions characterizing the physical properties 
(spatial or spectral) of fractals make them a useful testing ground for dimen- 
sionality dependent physical problems. This paper addresses specific problems 
including the behavior of the heat kernel and spectral zeta functions on frac- 
tals and their importance in the expression of spectral properties in quantum 
physics. Finally, we apply these results to specific problems such as thermo- 
dynamics of quantum radiation by a fractal blackbody. 



1. Introduction 

The interest in the behavior of fractals (a word coined by B. Mandelbrot in 
the 1970's [T] but without well agreed definition) goes back to the study by math- 
ematicians of strange objects hardly defined by their topology such as the Koch 
curve or the Sierpinski gasket. These objects are described by continuous but not 
differentiable functions. At about the same time, probabilists (Levy, Wiener, Doob, 
Ito, Kolmogorov) and physicists (Smoluchowski, Einstein, Perrin and Langevin to 
name a few) have been working to give a basis to the theory of brownian motion, 
yet another example of fractal object. More recently, physicists have recognized 
the ubiquitous character of fractals in almost all field of physics, including complex 
condensed matter ([2]- [8]), phase transitions [9lll0|. turbulence quantum field 
theory ( |12j - [l6] ) and aspects of stochastic processes [T] , 1171 [18] . An important ef- 
fort span over more than two decades led to new ideas and concepts to characterize 
fractals. Notions of self-similarity, iterative maps, fixed points and the identifica- 
tion of distinct fractal dimensions to characterize basic physical properties have 
been instrumental in the understanding of these objects. 

Since the early 1980's, mathematicians have opened new important directions 
by being able to define properly brownian motion on some classes of fractals ([19]- 
]) and Laplacian operators on these structures |24|, I23j . Progress along these 
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two directions have led to a vast literature and it would be a hopeless task to list 
it exhaustively. Most of these results have been summarized in textbooks {[24\- 
[26] ) and reviews [27 . These progresses prove to be instrumental in physics since 
they allow to go beyond phenomenological scaling relations towards a quantitative 
analysis of fractal structures. 

A useful description in field theory is provided by path integrals and more 
generally functional integrals 28j . Eventually, their evaluation, within one- loop 
approximation, boils down to the calculation of determinants of operators {e.g. 
the Laplacian) and their expression in terms of spectral functions among which the 
most useful and popular are the heat kernel and corresponding zeta functions (l30]- 
[32j ). Their evaluation on Euclidean manifolds reveal instrumental in quantum 
field theory [46], in statistical mechanics and the theory of phase transitions to 
name a few j29| . The new tools provided by mathematicians allow to extend these 
approaches to fractal structures f |33| - [35] ). 

It is the purpose of this paper to give an account of some of these new results. 
We shall present examples in the realm of wave and heat propagation on fractals, 
quantum mechanics, quantum field theory and statistical mechanics. 

Why is it interesting to study physical phenomena on fractal structures ? Frac- 
tals define a useful test ground for dimensionality dependent physical phenomena. 
Indeed, many physical phenomena reveal being critically dependent upon space 
dimensionality. Relevant examples include Anderson localization, Bose-Einstein 
condensation, onset of superfiuidity (Mermin-Wagner-Coleman theorem). On Eu- 
clidean manifolds, there is a single space dimensionality so that it is usually not 
possible to identify the meaning of a dependence upon dimensionality. On a frac- 
tal, as we shall see, there exist distinct dimensions which account for geometric, 
spectral or stochastic informations. A purpose of this paper is to discuss problems 
where this dependence plays a role. Fractals also provide a new playground for well 
designed new experimental setups |36j . An example is provided by spontaneous 
emission from atoms embedded in structures whose quantum vacuum has a fractal 
structure [37j . 

2. Discrete scaling symmetry - Self similarity - Definitions 

As opposed to Euclidean spaces characterized by translation invariance, self- 
similar (fractal) structures possess a dilatation symmetry of their physical proper- 
ties, each characterized by a specific fractal dimension. To illustrate our results, we 
shall consider throughout this paper simple examples of fractals such as the Siepin- 
ski gasket, families of diamond fractals represented on FiglHor Cantor sets, but the 
reader should keep in mind that our results apply to a broader class of self-similar 
deterministic fractals. Much less is known about other related but more complex 
systems including random fractals (e.g. critical percolation clusters), multifractals, 
T-graphs, trees or more generally objects which do not exhibit an exact decimation 
symmetry like in deterministic fractals considered here |38) . 

A fractal is an iterative structure. Let us give a simple example of dimensional 
characterization. To that purpose, consider a triadic Cantor set obtained iteratively 
by removing from a given initial interval of length Lq the middle third part. Define 
a uniform mass density on this interval so that the initial mass is Mq for a length 
Lq. After the first iteration, the mass is 2Mo while the length is Sip- After the n-th 
iteration, the mass is M„ = 2"Mo while the length becomes L„ = 3"io; so that the 
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Figure 1. First 2 iterations of the diamond fractals D4 2J 2 
and Dq,3. Their respective branching factors (defined in the text) 
are B = 1,2, 1. 



fractal dimension describing the mass density, the geometric HausdorfF dimension 
dh, is given by, 

, , lnAf„ In 2 

(2.1) dh = lim„_^oc-i — T~ = 

In Ln In 3 

An alternative way to obtain this result, is to start from the scaling relation 
between masses M{L) at different lengths L, namely M{L) = Mq{L) + ^M{3L) 
for the triadic Cantor set, where Mq{L) is the initial mass. More generally, we are 
interested in the solution of the equation 

(2.2) f{x) - g{x) + i/(ax) 

with determined scaling parameters a and h and a given initial function g{x). This 
form defines a discrete scaling symmetry as opposed to continuous scaling which 
extends (|2.2p to any rescaling a of the variable x. For the specific case 5 = 0, a 
solution of (|2.2p can be sought under the form f{x) — x°'F{x) for generally non 
constant F{x). Inserting into (|2.2p . and choosing a — ln6/lna, we obtain that F 
must fulfiU F{ax) = F{x). Redefining F{x) = G'(lnx/lna) leads to G(lna:/lna + 
1) = G(lna;/lna), namely the general solution of (12. 2p is of the form. 



(2.3) /(a;) =a;'" 



Inx 
In a 



where G is a periodic function of its argument of period unity. The scaling form 
can be iterated so that it rewrites 



(2.4) /(x) = ^ 6-"5(a"2:) 



ri=0 



This form makes explicit an important feature of discrete scaling symmetry namely 
the scaling function f{x) can be written as a series with exponentially growing 
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coefficients b " and a" instead of polynomial growth. This is an essential and 
ubiquitous feature of self-similarity that we shall encounter all along this paper. 

2.1. Mellin transform. To study functions involving dilatation scaling sym- 
metry, it is desirable to find an appropriate transform equivalent to Laplace or 
Fourier transforms for translation symmetry. The Mellin transform plays this role 
[39j . For a function f{t) defined on the positive real axis < t < cx), its Mellin 
transform 

/>oo 

(2.5) Mj{s)= / dtt'-^f{t) 

Jo 

is defined on the complex plane. The Mellin transform of f{t) is the Laplace 
transform of the function g{t) = f (e^*), a property which precisely accounts for the 
dilatation symmetry in (|2.2|) or the exponential behavior in ()2.4|) . An important 
property of the Mellin transform can be stated as follows. If f{t) is analytic in 
< t < oo and f{t) 0(t"") for i and f{t) = 0{t-^) for t oo, with 
a < P, then the Mellin transform Mf{s) is analytic in the strip a < s < /3, and 
fit) = (l/2i7r) Mf{s)t-' ds where a < a < p. Another suitable expression 

of the Mellin transform is provided by the zeta function of the function /(x) defined 

by 

1 

(2.6) Q^s)^—J^ dxx^-'f{x) 

which takes into account suitable properties of the Euler F function as we shall see 
later. The inverse transform is thus 

-1 />CJ" + ZOO 

(2.7) /(^)^ / dsx-^T{s)Cf{s) 

A direct calculation of the zeta function of a function f{x) with a discrete scaling 
symmetry (j2.2[) is 

(2.8) C/(.) = ^-f^ 

The behavior of f{x) is driven by the poles of (f. Disregarding at this stage the 
pole structure of Cg, we focus on the poles in (|2.8p which result from the scaling 
symmetry. There are the solutions, s„, of foa" = 1 namely, 

, , In (1/5) 2iTTn 

(2.9) s„ = + 

In a in a 

for integer n. The origin of these poles is a direct consequence of the exponential 
behavior of the coefficients in (12. 4p . These complex poles have been identified 
with complex valued dimensions of self-similar fractal systems |26) . By an inverse 
transform, we have 

^o^r^\ ft \ -l^iM \^ n n /■ / n ~2i^n('^)_ _ IsiiZid „ / 27rlni 

(2.10) f(X)=X 1»» F(s„)Cs(Sn)e Vlna^^j; l„a Q I ——— 

n— — oo ^ 

namely, the general solution of (|2.2|) is the product of a power law behavior char- 
acterized by a fractal dimension In (l/b) / Ina times a periodic function G{x -I- 1) = 
G{x) of its argument 27rlna;/lna. This log-periodic behavior is tightly related to 
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the existence of a discrete scaling symmetry as expressed by p.2|) and it consti- 
tutes its fingerprint together with the power law prefactor determining the fractal 
dimension. 

Log-periodic functions have a long history ^40j including the well known Weier- 
strass function determined by the series, 

oo 

(2.11) W^(t) = ^ a" cos (6" i) ■ 

They appear in a broad range of problems such as renormalization group and the 
theory of phase transitions (91 llOj , Markov processes in complex media [18 , , tur- 
bulence and fractals [6l [7] . Let us also mention that the scaling relation (12. 2p 
with definite parameters (a, b) is a particular example of more general iteration 
processes described by the Poincare equation f{az) = V {f{z)) where 'P{x) is a 
polynomial [40j . This is actively studied in the mathematical literature [27.. 



Remark 2.1. The periodic function G{x) becomes constant when all residues 
in the inverse Mellin transform vanish except for n — 0. This is the case for 
instance for the mass density M(L) of the interval of length L or generally of any 
Euclidean manifold of dimension d, where M{L) (x is expected without log- 
periodic behavior. This can be checked either by a direct calculation or by noticing 
that if the scaling symmetry M{aL) = bM{L) takes place for any value of a and 
not only for a fixed one, then, averaging over a washes out the oscillations leaving 
only the expected power law behavior. 

2.2. A variational derivation. It is of interest to present another derivation 
of the form (|2.10p not based on the inverse Mellin transform. This provides an- 
other point of view which is useful in cases where the pole structure cannot be easily 
retrieved. To that purpose, we start from the iterated form (|2.4|) of f{x) and to esti- 
mate it, we use a saddle point approximation to find the value of n which dominates 
f{x). By differentiating with respect to n, we obtain —{\nb)g + {\na)a^ x g' = 0. 
Defining u = x, we have — (ln6)g(u) -I- {In a)u x g' (u) = 0. This equation in the 
variable u admits a solution N such that u — x and then, 

,T luii — Inx 

(2.12) . 

In a 

N is not necessarily an integer, namely N — uq + t, where uq is an integer and 

< < < 1. Then, we write f{x) in (|2.4|) under the form, 

(2.13) 

OO OO 

m— — no + l m— — no + 1 

Using (I2.12P leads to b^^ = (u/a;)'"'"^^^''^'"", and since, u and a are fixed quantities, 
then no — )■ — oo for x — >■ 0, so that 

ln(l/i)) OO 

(2.14) /(•^) = (^) E fe-'™-*^ 5 



m— — oo 



where the series is now a periodic function of t of period unity, namely using ()2.12|) , 
a periodic function of In x/ In a of period unity. We thus recover ()2.10|) . 



6 



ERIC AKKERMANS 



3. Heat kernel and spectral functions - Generalities 

The study of spectral functions on manifolds has a long and very successful his- 
tory which traces back to the early XXth century with considerations put forward 
by Lorentz on the blackbody radiation |41j . Given an Euclidean manifold, it is 
possible to retrieve some of its geometric characteristics such as its volume, surface, 
curvature, Euler-Poincare characteristics by studying the spectral properties of the 
Laplace operator defined on the manifold [41j, \42\ I49j . In other words, the Lapla- 
cian or more precisely its eigenvalue spectrum can be viewed as a "ruler" which 
allows to span the manifold. Since stationary wave and heat equations are both 
governed by the Laplacian, looking at the heat flow or at (scalar) electromagnetic 
wave propagation are equivalent ways for extracting geometric characteristics of 
Euclidean manifolds. This has led to seminal works in mathematics starting with 
H. Weyl [43] and the expansion which bears his name culminating in the celebrated 
question of M. Kac, "Can you hear the shape of a drum ?" ^44; and the negative 
answer provided by examples of isospectral domains of different shapes |45j . 

The direct relation between spectral functions of Laplace operators and path 
integral in quantum field theory [46] or equivalent functional forms of the partition 
function in statistical mechanics are the underlying reasons for the enduring success 
of methods based on spectral functions of Laplace or corresponding Dirac operators 
in gauge theories |47j . 

Let us first illustrate these ideas using simple examples (see Chapter 5 in [50,). 
Consider first the diffusion of heat ip{x, t) along the infinite, unbounded real line 
(d = 1). The corresponding heat equation is 

(3.) '-±-D^, 

where the diffusion coefficient D sets units of length and time. The Green's solution 
of this equation is 

1 ijc-y)^ 

(3.2) Pd=iix,y,t) 



1/2' 



{4nDt) 

In a probabilistic interpretation, Pd=i{x^y,t) represents the probability density for 
a particle to diffuse from an initial position x to a final position ?/ in a time t. Its 
generalization Pd{x,y,t) to the c?-dimensional free space is obtained from (|3.2|) by 
replacing the exponent 1/2 in the denominator by d/2. A way to characterize the 
space geometry which sustains the heat flow, is to consider the heat kernel defined 

by 

(3.3) Za{t)^ [ d''xPd{x,x,t) = 

where the integral is over a volume V defined qualitatively without yet specifying 
boundary conditions. The heat kernel thus defined is an integral over all closed 
diffusing trajectories (starting and ending at a point x) within the volume V. The 
Green's function Pd is obtained for t > 0, from the normalized eigenfunctions (j)n{x) 
and non negative eigenvalues A„ (with degeneracy (?„) of the heat equation (j3.1[) in 
d dimensions, as 



(3.4) Pd(x,y,i) = ^5„0„(a;)0;(2/)e" 



A„ t 
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SO that 

(3.5) Z,(0-^5«e-^"* = Tr(e^*^) . 

n 

This relation between the heat kernel Zd{t) and the Laplace operator — A, expresses 
the distribution of closed diffusive trajectories within a manifold in terms of the 
spectrum of the Laplacian. This relation for the heat kernel is instrumental in 
calculating Euclidean path integrals |46| , partition functions |29j and other spectral 
and transport quantities |50j . 



3.1. The Weyl expansion. Expression (j3.5|) is also very convenient to eval- 
uate Zd{t) for manifolds with boundaries. Consider the simple case of diffusion 
on an interval of length L. The corresponding eigenvalue spectrum of the Lapla- 
cian (we set D = 1 for convenience) is given by A„ = (n-KjVf' where n is an 
non zero integer for Dirichlet {T>) boundary conditions, (\){x = 0,t) = (f>{L,t) — 0, 
whereas it includes the zero mode n — for Neumann (Af) boundary conditions 
d(f>{x,t)\x=o — d(f>{x,t)\x=L — 0. The corresponding heat kernels Zjsf,v{t) are thus 
related by 

oo 

(3.6) ZAA(i) = ^e-("^/^)'* = 1 + Zi,(t) 

The use of the Poisson formula allows to write the small time asymptotic expansion 

(3.7) Z^^T>{t) = ^ 



This is the simplest example of a Weyl expansion. For a two-dimensional domain of 
arbitrary shape with surface S and boundary length L, we have the Weyl expansion 
corresponding to Dirichlet boundary conditions [51| . 



where the constant term, 1/6, results from the integral of the local curvature of the 
boundary. More generally for a d-dimensional Euclidean manifold of hypervolume 
V, hypersurface S, etc., the Weyl asymptotic expansion (restoring the diffusion 
coefficient D), involves powers of 1/ (47r ZJi)'-'*"*''^^: 

(3-9) Zd(t) ~ -— ^— - ad ^ 



where ad is a constant which depends on boundary conditions |32|, \41\ I42j . The 
Weyl asymptotic formula provides a small time expansion for Zd{t). Physically, 
it describes the behavior of a diffusive particle initially released at some point in 
the manifold. At small time, it experiences a free space diffusion insensitive to the 
boundaries (volume term). At later times, the particle starts feeling the boundary 
(surface term), then its shape (local curvature term), etc.. 

To conclude this part, we discuss a point concerning the dispersion relation 
{i.e., the relation between time and length units) which will prove relevant when 
considering diffusion on fractals. Zd is, as defined in (|3.5p . a dimensionless quan- 
tity. On the other hand, the Laplacian has dimensions of the inverse of a length 
squared which allows to retrieve from the heat kernel geometric information about 
the manifold. We thus need to insert a diffusion coefficient D in order for V^^'^/D 
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in p.9|) to have units of time. The diffusion coefficient D expresses the underlying 
physics of the diffusion flow, and it is related to the relevant sources of diffusion 
by an appropriate phenomenological relation [e.g. temperature T and viscosity tr 
where D — RT/6n rj aj\f = ksTa) generally known as the Einstein relation 56]. 
For instance, in the specific case of a covariant diffusion equation as studied in quan- 
tum mesoscopic physics [50j , or in Euclidean time formulation of the Schrodinger 
equation of a particle of mass m, we have D = h/2m. 

3.2. Spectral determinant - Density of states and spectral zeta func- 
tion. There is an important and useful relation between the heat kernel, its ge- 
ometrical content and the density of states p{X) of the Laplacian defined on the 
corresponding manifold. To find it, we define the spectral determinant associated 
to the eigenvalue spectrum A„ of —A, 

(3.10) S{j)^det{-DA + j) = Y[{Xn+j) 

n 

where 7 is a real number. From the relation (13. 5p . it follows that 

(3.11) ^(^K-^^-E^=a;l^^^(7). 

The density of states p{X) ~ ^{^^ ^n) is thus directly related to the heat kernel 
and the spectral determinant through (see for instance Chapter 5 in [50j ) 

(3.12) p{\) = -ilim„^o+Ini-f-lnS'(7) 

TT d-f 

where 7 is now complex valued, 7 = —A + irj. This last relation proves useful but 
uneasy to implement since the spectral determinant £'(7) is defined by the product 
of eigenvalues A„. This product is infinite, so that its definition is formal. To give 
an interpretation to S'(7), we resort to the spectral Ca function, associated to the 
Laplacian and defined by 



(3.13) Ca(s)-Ei7 



n " 

This function is well-defined for all values of s for which the series converges. Using 
the identity 

(3-^4) A^ ^ i ' 

we write 

(3.15) Ca(s) = — / dtt^-iTr(e^^*) . 

-L Jo 

Ca(s) thus defined, is the Mellin transform (|2.5p of the heat kernel. It is convergent, 
and its analytic continuation in the complex plane defines a meromorphic function 
in s which is analytic at s = 0. We use this analyticity and the identity : 

(3.16) :^A-^|,=o = -lnA„ 

as 

to express the spectral determinant as 

(3.17) ln5(7 = 0) = -^CA(s)U=o 
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which is well defined. 

As an example, we consider the Laplacian on an interval of length L with 
Dirichlet boundary conditions (and set D = 1). Inserting the expression A„ = 
{mr/L f of the eigenvalues into (|3.13p . we obtain 

(3.18) Ca(.) = E(^) U2s). 

The Riemann zeta function Ci?(2s) has a simple pole at 2s = 1 and more generally 
a simple pole at 2s = d for a d-dimensional Euclidean manifold. Thus, the inverse 

Mellin transform (1^ provides directly Z-pit) = (L/27r) r(l/2)t"i/2 ^ namely 

p.7|) . Using p.l2p . we deduce the main contribution to the corresponding density 
of states 

(3.19) pd.i(A) ^ — ^ 

27r V A 

and from (j3.9p . we obtain the generalization of these results to d-dimensional Eu- 
clidean manifolds. 

To conclude this part, we note that the short time Weyl expansion of the heat kernel 
is related to the pole structure of the zeta function Ca associated to the Laplacian. 



3.3. Counting function - Spectral zeta function - Wavelet transform. 

We consider the counting function N{X) defined by 

(3.21) NiX)=J2l. 

A„<A 

The counting function A^(A) and the density of states p(A) — '^^J\^^ — X]n'^('^ ~ 
A„) are the sum of step 9 and Dirac S functions characterizing the distribution 
of eigenvalues A„ of the Laplacian. The density of states p(A) is related to the 
spectral determinant p. 121) . to the heat kernel and to the spectral zeta function. 
Unlike the counting function, the density of states is often ill-defined since A^(A) is 
not differentiable {e.g. for a fractal spectrum). 

Consider the number C(A) of eigenvalues of the Laplace operator —A within 
the interval [0, A] included into the support J of the whole spectrum. It is given in 
terms of the counting function as 

(3.22) C(A) iV(A) - iV(0) = J dN{e)e{X-e) 

where the integral extends over the whole spectrum J with the spectral measure 
dN{e). By definition p.5p . the heat kernel Z(t) can be expressed as (with D = 1) 



(3.23) 



Z{t) = Tr (e* = J dN{X) e'* ^ 
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The spectral zeta function Ca(s) defined in p.isp thus becomes, 
(3.24) = j dN{\)\-' 

as a result of the identity T{s) = dxx^~^ e~^. It is thus possible to relate the 
Mellin transform Mc{s) of C(A) to Ca(s), namely, 



/ dC{X)X-' = f X^'df f 
J J J J \Jj 



Mc{s) = / dC{X)\-'= / \-'d[ / dN{e)e{\-t) 

A"" / dN{t)6{\~e) 
I J J J 

(3.25) = Ca(s) 

so that the spectral function Ca(s) is the Mellin transform of the counting function. 
This relation is a particular example of a more general property of the Mellin 
transform which can be viewed as a convolution theorem. To state it, consider an 
integrable function g{x), suitably normalized to unity, Jg digit) — 1, and whose 
Mellin transform Mg{s) is defined. The quantity 

(3.26) Wg (Afc, t) = ^ dN{X) g{t\X- Ab|) , 

defined for i > 0, is called the wavelet transform of the counting function N{\). In- 
tuitively, the wavelet transform can be viewed as a mathematical microscope which 
probes the counting function at a point Ab with an optics specified by the choice of 
the specific wavelet g{x). An important property of the wavelet transform is that 
it preserves a discrete scaling symmetry ()2.2[) of the probed function. Performing 
a Mellin transform of Wg (Ab, t) w.r.t. the variable t gives 

(3.27) M [Wg (Afc, t)] = Ca (s, -Afc) Afg(s) 

where Ca (s, — Af,) = Jj dN{X) |A — Ah|^* is a shifted zeta function. From the last 
relation, it is immediate to obtain the relation between different wavelet transforms 
respectively specified by the "optics" f{t) and g{t) of Mellin transforms Mj and 
Mg, namely 

.o.n^ M[Wg{X,,t)] M[Wf{Xk,t)] ^ . , . 

(3-28) rj-r^ = l^rT~\ =CA(s,-Afc). 

Mg{s) Mf{s) 

An interesting application of this convolution rule is obtained by taking A;, = 
and the wavelet g{t) = such that Mg(s) = (2/7r)r(s - 1) sin7r(s - l)/2 for 

< Rc(s) < 2. We obtain\-om (jX^ . 

(3.29) M 



TT 



J 



,,r/, N sinXt 
dNiX) ^ 



-r(s-l) sin5(s-l)CA(s) 

TT Z 



SO that from the inverse Mellin transform, 



7r(s - 1) 



r(s)CA(s) 



Therefore different probes of the spectrum of the Laplacian can be related one to 
another and expressed in terms of the spectral zeta function. This is a powerful 
result with interesting physical consequences since distinct physical phenomena 
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usually involve their own probing function {e.g. the Fermi golden rule describing 
the response of the spectrum to an external perturbation (for instance spontaneous 
emission), involves the sine probe g{t)) |37j and all are related to the heat kernel 
characterized by an exponential probe. 

4. Laplacian on fractals - Heat kernel and spectral zeta function 

This section is devoted to studying the heat kernel on fractals, namely on 
systems whose geometry is characterized by a discrete scaling symmetry as defined 
in (|2.2[) . In other words, we would like to extend the previous analysis and the 
relation between geometric and spectral properties to fractals. This is a vast subject 
and we do not intend to be exhaustive but rather to study specific but generic 
enough examples in order to highlight some of the more salient results and open 
questions. 

We have seen that geometric characteristics of Euclidean manifolds can be re- 
trieved from the spectrum of the corresponding Laplace operator — A. The relevant 
spectral tools are the heat kernel (|3.5[) and the spectral zeta function Ca(s) defined 
by (|3.15p . How do they generalize to fractals ? Geometric information about a 
fractal is characterized by its Hausdorff dimension dh- On the other hand, the very 
notion of volume or surface of a fractal is rather ill-defined. Then, is it possible to 
define a heat kernel and a Weyl expansion for fractals and if it is so, what kind of 
information does it provide. Those questions have been and still are in the focus 
of intensive works from the mathematical community where an abundant number 
of important results has been proved ([19j-[27]). It is not our purpose to review 
them, but rather to show their relevance and usefulness in physics. An important 
step has been to prove, using either a probabilistic or an analytic approach, that 
a Laplace operator can be properly defined on fractal structures [24j as well as 
corresponding heat kernel and spectral zeta-function |20| . 

As we have seen in Section 2, the characteristic feature of a discrete scaling 
symmetry is the existence of a tower of complex poles in the zeta function associated 
to the relevant quantity. These complex poles are a direct expression of the scaling 
form (|2.2p and of the exponential behavior of the coefficients in the iteration series 
(|2.4|) . Are there similar properties of the spectral zeta function of the Laplacian on 
a fractal. The answer to this question is positive. There exist two characteristic 
parameters (a, h) = {l'^™ , Z'^''), such that the heat kernel of the Laplacian on a self- 
similar fractal obeys the scaling relation ()2.2p . These scaling parameters involve 
the geometric Hausdorff dimension dh, the walk dimension d^ that we shall define 
and discuss later on and an elementary step length / describing size scaling. 





dh 


dw 


ds = '2dh/d^, 
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Figure 2. Fractal dimensions and size scaling factor for diamond 
fractals and for the Sierpinski gasket. For 1)6,2, the spectral di- 
mension is ds ~ 2.58, and for -De, 3, ds ~ 1.63. 
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To proceed further, we consider the specific case of a diamond fractal (0) 
(see Figs H] and HI). It has been shown 33] that the corresponding heat kernel 
Z(y{t) has a closed expression of the form (|2.4I) where the initial function to be 
iterated is the heat kernel Z-pit) of the interval of unit length with Dirichlet (I?) 
boundary conditions, given by dSH), namely, Zv{t) = XlfeLi e"'' * (we set £> = 1 
for convenience), so that, 

oo 

(4.1) Z(, ^Zv{t)+Bj2 Lt Zv {Lf- t) 

71=0 

where we have defined the total length i„ = I" of the diamond fractal at step 
n of the iteration. The coefficient B = (^l'^'^^^ — l) is the branching factor of 
the fractal (see FiglTJ and the integer Z'*'* is the number of links into which a 
given link is divided. As in (|2.4I) . the series for Z(y{t) involves coefficients which 
behave exponentially with the iteration n. We therefore expect an expression of 
the form (|2.10p with log-periodic oscillations with time. To show it, we calculate 
the corresponding spectral zeta function, 

(4.2) Co = T^/ dtt^~'Z^{t). 
An elementary calculation |33l leads to 



d^jjS 



n=0 



(4.3) ^ C«^^.._,n-^ 



1 _ [dh-d^^s 



where Cr{2s) is the Riemann zeta function. Note that a similar structure exists for 
the Sierpinski gasket [20|, I23j , with the Riemann zeta function factor replaced by 
another zeta function. Coi^) complex poles given by 

dh ^ 2i7rm _ ds ^ 2i7rm 
dw c?u, In Z 2 In Z ' 

where m is an integer. These complex poles have been identified with complex 
dimensions for fractals j26] . The fractal dimension dg = 2dfi/dw is called spectral 
dimension. It has been obtained in earlier works in the physics community [4l |5] 
and recognized as the relevant fractal dimension (unlike dh) underlying spectral 
properties of self-similar fractals. 

From the inverse Mellin transform (12.71). we obtain 



V 

^OW = ^ E r(s™)Ci?(2s™)A2^"e2'^""^*/('^™'"') 



m— - 



= ^^Hrf^Th^ 

where Go is a periodic function of its argument of period unity. These log-periodic 
oscillations are represented on FiglS] and we note that the higher order (with m) 
complex poles give much smaller contributions, a result related to the steep decrease 
of the Euler F function. A similar behavior has been found numerically for the 
Sierpinksi gasket 
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Figure 3. Heat kernel Zoit) at small time, normalized by the 
leading non-oscillating term, for the fractal diamond 1)4,2 ■ The 
solid [blue] curve is exact; the dashed [red] curve is the approxi- 
mate expression (|4.5p . At very small t, these curves are indistin- 
guishable, as shown in the inset plot. The relative amplitude of 
the oscillations remains constant as t — ?> .33 . 



For mathematical discussions of oscillations in heat kernel estimates see |52l 
I53j . In particular, there is a class of fractals where oscillations are related to large 
gaps in the spectrum. This topic is a subject of active current research |54|, I55[ 
and references therein] . 

4.1. Spectral volume - Anomalous walk dimension. We have seen that 
for a d-dimensions Euclidean manifold, the dominant contribution to the Weyl 
expansion takes the form p.9|) . i.e., it is proportional to the volume V of the 
manifold. This dependence, as has been emphasized, is a consequence of the fact 
that the Laplacian sets units of length. In other words, the quantity i'^/^Tr (e~'^) 
has a well defined limit, V, for t — ^ 0. This limit can be viewed as ^ = L** where L 
is the characteristic length set by the Laplacian. There is another consequence of 
the form L'^ jt'^l'^ of the dominant term of the Weyl expansion p.9p . Since Z{fy is 
dimensionless, this imposes the usual Euclidean dispersion (xt characteristic of 
a diffusion processes where units are matched by means of the diffusion coefficient 
D. This can be checked directly from the differential equation p.ip or from its 
solution p.2[) . But for diffusion on a fractal there is no similar local equation nor 
closed expression of the local propagator P{r,r',t). It is nevertheless possible to 
express an equivalent dispersion which is solely a property of the Laplace operator. 

To that purpose, we note from (|4.5p . that the quantity f^'^^ Zq/Gq has a finite 
limit for t — > which we denoted Vg. Restoring length and time units (through 
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D) and comparing to (|3.3p for the Euclidean case, leads necessarily to the form 
Vs = if'' where Ls is the characteristic spectral length associated to the Laplacian 
on the diamond fractal and is the Hausdorff geometric dimension. Inserting this 
expression into Z(-, allows to rewrite the power law leading prefactor under the form 



2^2(ih/d, . Since it is a dimensionless quantity, this implies that the disper- 

sion for fractals is of the form oc This new dispersion plays an essential 

role in the characterization of diffusion on fractals. It is a constitutive equation 
usually written under the form of a characteristic mean square displacement law, 

(4.6) (r2(t)) oc t^/'*™ 

where the exponent dw = 2dh/ds is called "anomalous walk dimension". Thus, the 
expression, 

for the heat kernel on a fractal takes a form analogous to its Euclidean counterpart 
p.ip but where the geometric volume V = L'^ is now replaced by a spectral volume 
Vs = if'* which results directly from the spectral properties of the Laplacian. It 
is to be noted that the geometric volume of a fractal becomes infinite so that it is 
unlikely to appear in the Weyl expansion of Z^/Go which has a finite limit 
for t 0. 



5. Thermodynamics on photons : The fractal blackbody |34| 

5.1. Thermodynamics of the quantum radiation - GeneraHties. A di- 
rect application of the considerations developed in the previous sections is to the 
study of statistical mechanics of quantum radiation (photons) in fractal structures. 
Indeed, a basic aspect of it, namely the blackbody radiation precisely addresses 
the relation between the electromagnetic modes inside a cavity (a manifold) and 
its geometric characteristics. Historically, the seminal work of H. Weyl which led, 
among other results, to the expansion p.9p has been motivated by considerations 
raised by H. Lorentz about the dependence of the Jean's radiation law upon the 
volume of the cavity at the exclusion of other geometric characteristics. 

Although the derivation of the thermodynamic partition function of quantum 
radiation constitutes well known textbook materials, we wish to re-examine it to- 
wards its application to fractals in order to emphasize some key points and basic 
assumptions. The purpose of this derivation is to show that the partition function 
is directly related to the spectral zeta function so that it can be calculated for Eu- 
clidean and fractal manifolds as well. Moreover the geometric information about a 
manifold retrieved from thermodynamics relates directly to the spectral geometry 
of the Laplace operator. This point is not always emphasized in textbooks which 
prefer to consider combinatorics of mode counting in simple cubic geometries, an 
approach which relies heavily on the existence of Fourier transform and phase space 
quantization cells, a tool which is not available for fractals. 

For convenience, we consider thermodynamics of a scalar massless field (scalar 
photons) which simplifies the calculation due to the vanishing chemical potential 
resulting from zero mass. Generalization to massive fields is possible though slightly 
more cumbersome [59^ 160] . The spectral partition function (including the zero 
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mode Casimir contribution) of a quantum oscillator of frequency u at temperature 
T = l/ksP is \nZ{T,uj) = - In (l - e'^'^). Using the identity 



(5-1) -H^-^-n=t[%e--'^h 

together with the Poisson formula 

+00 I — +00 



2 2 / 
g-7r n It 



leads to 

(5.3) 



1 /'OO J +00 



The operator 9q = df defined with periodic boundary conditions cf>{t + hf3) = (pit) 
admits a discrete spectrum M known as Matsubara frequencies 27rn//i/3, so that 

+00 

(5.4) ^ e-(^)'^ = TVm (e"^^°) • 



n— — 00 

Identifiying = c^k"^ with the eigenvalues of c^A, and tracing over all modes, 
we obtain for the partition function of the quantum radiation at temperature T = 
l/kBp, 

(5.5) In Z{T, y) = \ ^Tr^ (e""^'^) Ttm (e"^' 

where Tvm denotes the trace over the modes of the Laplacian defined on the mani- 
fold Ai which contains the anticipated dependence on the volume V . Finally using 
the identity InC = - ^e"'^'", and InDet O = Tr InC, we obtain the elegant 
and compact dimensionless form for the partition function, 

(5.6) \nZ{T,V) = -ilnDetMxM (4' + ^' a) 

(5.7) = ^~^I;Mlv^(e-SA) 

where we have defined the photon thermal wavelength = hPc and the dimen- 
sionless operator dl = {H/Sf dl so that /(r) = Tim (e-^^o) = ^^^-00^"^^''"''^ 
is a dimensionless Jacobi 6'3-function. 

From this last expression of the partition function, we see that it is indeed 
determined by the heat kernel of the Laplace operator on the manifold M. 

(5.8) /Ca<(t) =TrA<e-^^^^. 

This form is useful to describe in a closed way the thermodynamics of quantum 
radiation in a large hypercube of volume L'^, in d space dimensions. There the 
modes k are quantized in units of 2Trn/V^^'^ where n is a vector with integers 
components defining elementary cells in the reciprocal phase space. The previous 
expression of the partition function becomes 

(5.9) lnZ{0,V = L'^) = \nZ{hpcV-^l'^) = -^InDetMxM (^^ + L^F-^/'^a) , 
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where —A is a dimensionless Laplacian. Therefore Z is a function of the single 
variable LpV~^^'''. Standard thermodynamic quantities fohow immediately from 
this specific scahng behavior. For example, the equation of state, PV = U/d, 
relating the internal energy U of the radiation to its pressure P and the volume V, 
follows immediately from 

(5.10) U 
and 

(5.11) P 

The Stefan-Boltzmann law for the internal energy J7 is a consequence of the equation 
of state and the thermodynamic relation (f^)^ = T (^)^— P, while noticing from 
([5TT1) that P depends on T only, in the thermodynamic limit. We then obtain 
U = aVT'^'^^ where a is a constant to be determined from (|5.9|) . It is already 
apparent from this simple case that the volume dependence in the thermodynamic 
equation of state, PV = U/d, comes from spectrum of the Laplacian. 

For black-body radiation associated to Euclidean manifolds of complicated 
shape, it is difficult to make an explicit mode decomposition and find an ex- 
plicit expression like ()5.9p for the heat kernel. However, we can learn about the 
thermodynamic [large volume] limit from the Weyl expansion p.9|) of the heat 
kernel. Note that the large volume limit corresponds to V ^ LJ^, which is a 
"high temperature" limit ksT ^ hc/V^^'''. Keeping only the dominant volume 
term in (j3.9p . expression (|5.7p leads immediately to the familiar thermodynamic 
expressions [56] previously derived: InZ = (y/L^)Ci?(d l)r (^) /tt^'^+i^/^, 
P = {kBT/LfjCR{d + l)r (^) /7r(^+i)/2. Away from the thermodynamic limit, 
subdominant terms in p.9p lead to corrections that depend on the exact geometry 
of the volume enclosing the radiation, but the equation of state PV = U/d is always 
valid |42| . This formulation in terms of heat kernel makes it clear that the Weyl 
expansion is directly related to the thermodynamic limit of a black-body radiation 
system, so we can use the leading term as a definition of the volume probed by the 
photons as they attain thermal equilibrium. 

5.2. Thermodynamics of the quantum radiation on fractals. Based on 
previous results and especially expression (j5.7p of the partition function, we are now 
in position to study thermodynamics of quantum radiation on fractals [34]. The 
heat kernel /Cj^(t) equivalent of (|5.8I) but on a fractal T is obtained using (14. 7p . so 
that 

The thermodynamic partition function In is a function of the single 

variable L^T4~^^'^'', where v}^"^*" — Lg is the spectral length associated to the 
Laplacian on the corresponding fractal J-. Thermodynamic quantities and relations 
follow immediately from this specific scaling behavior and from (|5.10p and (|5.1ip to 
give the thermodynamic equation of state on a fractal as PVs = U /dh- The second 
important conclusion is that the actual expressions for pressure P, internal energy 
U, etc... will be modified on a fractal, not only by the appearance of the spectral 



5/3 



lnZ(r, V)^- 



1 /91nZ 



dV 



d\nZ{x)\ hcV-^/'^ 



dx 



Vd 
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dimension ds and spectral volume Vs, but also by the appearance of oscillatory 
terms arising from the behavior of the log-periodic function Gjf in ()5.12p . 

5.3. Vacuum Casimir energy. Another straightforward consequence of ()5.12|) 
is the expression of the zero temperature free energy namely the vacuum Casimir 
energy calculated in various complex geometries and recently on quantum graphs 
[57j . A general expression for the Casimir energy is obtained from the inverse 
Mellin transform of the spectral zeta function which gives 

InZiT.V) . -j(f)c.(4) 

(5.13) +1, f r(2s)Cfl(2s + l)CA^(s)ds 

The first term gives the standard zero temperature "vacuum Casimir energy" con- 
tribution f58|, proportional to (^m (~^)- On a fractal we have L = Lg and the 
previous expression leads for the vacuum Casimir energy to: 

Ihc f l\ 

^'■''^ ^° = 2i:^-b)- 

6. Conclusion and some open questions 

We have presented general features of deterministic self similar fractals which 
possess an exact decimation symmetry. To describe quantitatively diffusion pro- 
cesses and wave propagation on these structures, we have defined and develop 
spectral tools related to the Laplacian, such as the heat kernel and the spectral 
zeta function. After reviewing some of their main features for Euclidean manifolds, 
we have calculated them on fractals. This has enabled us to single out a number of 
specific features of fractals such as the spectral dimension dg, the spectral volume 
Vs, the existence of log-periodic oscillations of spectral quantities, an unusual dis- 
persion characterized by the walk dimension d^ = 2dh/ds, dh being the geometric 
Hausdorff dimension. These features show up when extending results of statistical 
mechanics and quantum field theory on fractal structures. To see this at work, we 
have considered the problem of thermodynamics of quantum radiation on fractals 
in the simplified scalar version as well as the vacuum Casimir energy. The results 
we have obtained generalize straightforwardly to related problems such as massive 
bosons. In that case, it is a direct consequence of the form of the heat kernel that 
Bose-Einstein condensation occurs for dg > 2 i.e. independently of the geometric 
dimension dh [59\ I60j . a fact already recognized long ago [5]. More involved is 
the problem of the occurrence of superfluidity (a phenomenon distinct from Bose- 
Einstein condensation) which depends on terms of higher order than the volume 
contribution in the Weyl expansion [60] . Beyond equilibrium situations, the prob- 
lem of quantum emission of radiation (either spontaneous or stimulated) is highly 
interesting either from a fundamental point of view or for applications. It can be 
shown [37] that the probability of emission, also known as vacuum persistence in 
the context of quantum field theory, is driven by another spectral kernel, the sine 
kernel (|3.30l) related to the heat kernel and which exhibits an analogous qualitative 
behavior. 

We have mostly discussed spectral quantities and not the behavior of non- 
diagonal terms in the propagator P{x,y,t) associated to diffusion. Very little is 
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known about it except for useful bounds 19 and a conjectured form supported by 
numerical results |61j . This part of the characterization of fractals is nonetheless 
essential to better understand stochastic processes, out of equilibrium statistical 
mechanics and large deviation physics. Preliminary results based on specific distri- 
butions of traps |18| or on an additivity principle [62] seem promising. 

Let us mention finally that the absence of translational invariance in fractals 
prevents using Fourier transforms. It is preferable instead to use Mellin trans- 
form as explained in section 2. This has important consequences in formulating 
an uncertainty principle relating direct to reciprocal spaces. We have emphasized 
that while the direct space is driven by the Hausdorff dimension dh, candidates 
for reciprocal spaces obtained from the Laplacian are instead driven by dg. This 
issue has far reaching consequences if one wishes to properly formulate canonical 
field quantization on fractal structures including spin. This leads to yet another 
challenge on fractals. A proper understanding of spin on Euclidean manifolds has 
required tools such as heat kernel and Weyl expansion of Dirac operators. Topolog- 
ical properties related to spin thus appear under the form of powerful results known 
as Index theorems [47j . Their extension to fractals requires a proper understanding 
of cohomology, connexions and definition of a Laplacian using Hodge theorem [63j . 

Acknowledgments. : It is a pleasure to thank G. Dunne, A. Teplyaev, E. 
Gurevich, D. Gittelman and O. Spielberg for collaboration and discussions. 
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